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ABSTRACT 


Basic  Firmer  models  are  tested  in  both  horizontal  and  vertical  wind 
tunnels  in  order  to  develop  dynamic  testing  techniques  and  to  perfect  new 
methods  of  non-linear  analysis.  The  simple  free  motions  of  pure  pitching 
and  pure  rolling  are  first  studied.  Then,  free  angular  motions  of  combined 
pitching,  yawing,  and  rolling  are  investigated.  While  the  linear  analysis 
methods  appear  to  yield  good  results,  the  new  non-linear  analysis  methods 
are  found  to  be  essential  in  providing  an  accurate  representation  of  the 
various  motions  and  in  providing  correct  values  for  the  various  static  and 
dynamic  stability  coefficients . 
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IN  i  ROuuuiluN 

The  Basic  Finner  was  designed  as  a  fundamental  research  configuration 
for  u-rguided  finned  missiles  following  World  War  II.  *  >2  Numerous  studies 
have  neen  carried  out  in  air  and  in  water.  The  air  tests  have  utilized  both 
wind  tunnels  and  aeroballistic  ranges  over  the  entire  speed  range,  subsonic, 
transonic,  supersonic,  and  hypersonic.  In  water,  tests  have  utilized  towing 
tanks,  rotating  arms,  water  tunnels,  and  free  drops.  Unfortunately,  the 
static  and  dynamic  stability  coefficients  obtained  from  the  various  testing 
techniques  and  facilities  do  not  always  agree.  Serious  unexplained  differences 
persist. 

There  is  a  fundamental  requirement  to  obtain  correct  values  for  the 
various  stability  coefficients  in  order  to  understand  and  predict  actual  missile 
performance  and  accuracy.  Therefore  the  introduction  of  a  basic  research 
program  to  develop  wind  tunnel  dynamic  testing  techniques  and  to  perfect 
methods  of  non  -  linear  data  analysis  was  initiated.  ^ 

The  Basic  Finner  is  first  tested  by  allowing  free  yawing  motion  about 
a  side  support  in  a  horizontal  wind  tunnel.  The  resulting  angular  position 
data  is  reduced  both  by  a  linear  method  of  data  analysis^  and  by  a  new  method 
of  non-linear  analysis. 

A  unique  vertical  wind  tunnel  is  used  in  the  studies  of  the  free  rolling 
motion  of  the  Basic  Finner  mounted  on  a  rear  sting  support.  Again,  both 
linear  and  non-linear  methods  of  analysis  are  employed.  In  these  studies 
particular  attention  is  directed  towards  the  determination  of  the  non-linear 
induced  roll  moment  coefficient (a)  ,  which  can  cause  Roll  Lock-In  and 
which  is  a  requirement  for  Catastrophic  Yaw. 

Special  emphasis  is  placed  on  testing  techniques  which  permit  complete 
model  motions  in  all  three  degrees  of  angular  freedom.  In  one  testing 
technique  the  Basic  Finner  is  mounted  on  a  side  support  in  a  horizontal  wind 
tunnel.  The  model  is  free  to  pitch,  yaw,  and  roll,  except  for  a  small  central 
section  which  does  not  roll.  In  another  3-D  testing  technique  the  model  is 
mounted  on  a  rear  sting  in  a  vertical  downdraft  wind  tunnel. 

In  all  cases  the  non-linear  methods  of  analysis  are  first  studied  by  us¬ 
ing  data  generated  on  a  high  speed  computer  by  numerical  integration  of  the 
exact  equations  of  motion.  The  results  obtained  by  applying  the  non-linear 
analysis  methods  to  this  computer  data  are  compared  with  the  known  com¬ 
puter  inputs.  The  non-linear  analysis  methods  are  then  applied  to  various 
dynamic  wind  tunnel  data. 
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In  these  studies  particular  attention  is  given  to  non-linearities  in  the 
static  moment  coefficient, Civ/jn  the  damping  and  lag  moment  coefficients, 
C  4 C  and  the  Magnus  moment  coefficient,  C.. 

Mpa 
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GENERAL  METHOD  OF  ANALYSIS 


In  general,  the  various  free  motions  of  a  missile  are  described  by  an 
equation  of  the  form 


off)  =Ke(X  +  iaj)t+a. 


It  is  desired  to  '‘fit”  this  equation  to  experimental  data.  The  above  equation, 
however,  is  non-linear  in  the  unknowns  CO  and  \  ;  therefore  the  convention 
al  method  of  least  squares  cannot  be  used. 

The  left  side  of  Equation  (1)  may  be  expressed  in  terms  of  the  experi¬ 
mental  data  asd(T)  +  V  >  where  v  is  the  residual  (the  difference  between 
the  calculated  point  and  the  data  point).  Equation  (1)  now  becomes, 


aCO+V  =  fttaT,K,X,w,S) 

The  five  unknowns  may  be  expressed  as , 

“t  =  aT„  +  AaT 

K  =  K„  +AK 

X  —  V  +AX 

co  —  ujo  +  Aw 

8  —  8n  +aS 


where  a T  >  K0  >  X0  »  etc.  are  initial  approximations.  Equation  2  now 
becomes  ’ 


a(t)W  =  fCtaTo+^aT,K0-t-AK,X0+AX,a!>+-Aw,8 +a8)<3> 


Expansion  of  Equation  (3)  in  a  Taylor  Series  gives, 


a(t)  +V  =  -l-AaT-||-  +AK  -|i- 

•4-AX  ^  +Aw  ^  +AS  jfe 


(4) 

If  th°  initial  approximations  are  good,  the  /\  terms  will  be  small  and 
the  higher  order  terms  may  be  neglected.  Equation  (4)  is  then  linear  in  the 
differential  corrections/^  a  T  ,/\  K  »  *  •  •  and  least  squares  methods 
can  be  applied  to  their  determination.  Once  the  differential  corrections  are 
determined  they  may  be  added  to  the  initial  estimates  of  the  unknowns  ciT  » 
K0  *  X0  »  •  •  •  and  the  process  may  be  repeated  to  determine  a  new  set 
of  differential  corrections.  This  Method  of  Differential  Corrections  is  re¬ 
peated  until  the  sum  of  the  squares  of  the  residuals  is  a  minimum  or  until 
there  is  no  significant  change  in  the  unknowns. 

The  Wobble  computer  program4  employs  the  Method  of  Differential 
Corrections  (MDC)  as  described  above  to  fit  the  solutions  of  the  governing 
differential  equations  of  motion  to  experimental  angular  data.  This  program 
is  equipped  to  fit  the  simple  equation,  Equation  (1),  or  the  more  complex 
solution  for  the  Quadricyclic  Theory 


(5) 


By  overlapping  fits  of  small  sections  of  data  the  motion  parameters 
can  be  determined  as  continuous  functions  of  time .  The  various  aerodynamic 
stability  coefficients  may  then  be  determined  from  these  parameters. 


SECTION  HI 


ONE  DEGREE  OF  FREEDOM  PITCHING  MOTION 


A  Basic  Finner  model  was  mounted  on  a  side  support  in  a  horizontal 
wind  tunnel  and  allowed  to  freely  yaw.  A  schematic  of  the  Basic  Finner  is 
shown  in  Figure  1. 


Figure  1.  Schematic  of  Basic  Finner 


In  one  case  ball  bearings  were  used.  In  a  more  optimum  case,  the  model 
was  supported  from  below  by  a  needle  resting  in  a  sapphire  jewel  cup.  These 
tests  were  all  conducted  at  low  subsonic  speeds. 

Linear  Analysis 

The  differential  equation  for  one  degree  of  freedom  pitching  motion  may 
be  written  as  3 


or 


Maa+M,<S+IV^+Mg£=I<fc  , 


a+N,  d  +  N2a  =  Nj  (7) 

The  general  solution,  when  the  coefficients  are  constants,  is  given  by 

a=  K^V  +  K3  (8) 

since  u),  and  \  =  \ 

t  *  i 


,  this  may  be  simplified  to 


where 


a-  a,-  +  Ke^cos(wT  +8) 


X,,  +  iwi?  -  ±i luF- 


r  ^  21u£ 

^  ~  “pTTsd 

CM?+CMd  ^  §JSD 


The  purpose  of  this  analysis  was  to  measure  precisely  the  angle  of 
attack,  a  ,  as  a  function  of  time  and  to  use  Equation  (9)  to  obtain  u)  and 


and  then  use  Equations  (11)  and  (12)  to  obtain 


and  v  Mn  1  v-w; 


It  was  also  desired  to  test  the  transformation  equations  for  the  stability 
coefficients  using  six  different  center  of  gravity  positions.  Once  the  moment 
stability  coefficients,  Cm„  and(CM.,+  Cm*)  have  been  measured  at  one 
center  of  gravity,  the  axis  of  rotation  may  be  changed  and  new  values  obtain¬ 
ed  for  each  of  the  six  positions  on  the  Basic  Finner  model. 

The  equations  which  relate  the  stability  coefficients  at  one  center  of 
gravity  (or  rotation  axis)  to  another  center  of  gravity  are^ 


r* 


:Ma0+  CZa 


“  2®  CMa  ‘  CM,+@C2 


The  sign  convention  for  the  forces  and  moments  discussed  above  are 
shown  in  Figure  2. 

A  model  of  the  Basic  Finner  was  mounted  on  precision  bearing  on  a 
side  support. 5.6  Optical  techniques  were  employed  to  obtain  the  angular 
data  and  measure  the  photographic  plates.  A  sample  of  the  1-D  oscillatory 
motion  is  presented  in  Figure  3. 
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The  results  of  the  UNIVAC  1 107  runs  in  applying  the  MDC  to  the 
experimental  uata  snuwcu  uic  excellence  of  the  data.  This  was  indicated  by 
a  small  probable  error  of  "fit":  the  average  was  .  14°.  This  is  an  independent 
indication  of  the  very  small  error  in  the  angular  measurement  of  the  model 


motion.  Also  the  values  of  the  stability  coefficients  Cm*  and(CMt|^CMA  ) 
were  well  determined.  The  average  accuracy  is  1 %  for  C  lrt  and  2%  for 


Figure  2.  Force  &  Moment  Diagram 


Figure  3.  Angle  of  Attack  versus  Time 
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These  moment  stability  derivatives  are  plotted  in  Figures  4  and  5  as 
a  function  of  center  of  gravity  position. 

The  static  moment  transformation  equation,  Equation  (13),  was  fitted 
by  least  squares  to  the  data  in  Figure  4.  Also  the  damping  moment  trans¬ 
formation  equation,  Equation  (14),  was  fitted  by  least  squares  to  the  data  of 
Figure  5. 

The  P.E.  of  "fit"  of  the  two  transformation  equations  is  of  the  same 
order  as  the  accuracy' of  Cm*  and  ( )  in  Table  I.  Thus  the  trans¬ 
formation  relations  were  tested  and  found  accurate.  We  have  also  obtained 
from  the  transformations  the  values  for  the  force  stability  derivatives  as  in 
Table  I. 
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Figure  4.  C{yj  versus  C.G.  Location 
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A  A  »*-•  A_#  A  •  A  »»••»  1W  4  •*— '  *»»»»*  »  A  *  V/  I  »  <  ■  A  «  « 

Static  Moment  Transformation  (C.G.  =5.0) 

Coefficient  Value  from  Fit  Probable  Error 

U*  -40.1  +.168 

C,  -12.9  +.126 

Probable  Error  of  Fit  +0.561 

Dynamic  Moment  Transformation  Data 

-439.9  +  12.5 

<-*,  *C„  -120.3  ±9.4 

Probable  Error  of  Fit  +  41 . 75 


By  applying  MDC  repeatedly  to  overlapping  short  sections  of  data,  non- 
linearities  in  the  stability  coefficients  may  be  studied.  In  Figures  6  and  7, 
the  stability  coefficients  are  plotted  as  functions  of  RMS  cs  and  RMS  q 
It  is  seen  that  Cm*  and  (  CMcj  +  ),  and  probably  and 

(  +  C*,;<)  are  quite  non-linear. 


RM.&  «  (<**$.) 


Fig.  6.  CM  ^  versus  RMS  «-< 


Fig.  7.  CMq  +  Cm*  versus  RMSq 

Application  of  the  linear  theory  to  overlapping  small  sections  of  data  gives 
only  an  indication  of  non-linearities  in  the  aerodynamic  stability  coefficients. 
Therefore  to  determine  more  accurately  these  variations,  a  non-linear 
mathematical  model  must  be  utilized. 

Non-Linear  Analysis  ^ 

For  the  case  where  the  stability  coefficients  are  nonlinear  functions  of 


angle  of  attack  as  follows: 

CM^a)  =  Cm^+  Cm^O?  -h . . 

C^)  -  CM„  . . . 


(15) 


'MOo  ~mCL 2 

The  approximate  solution  of  the  differential  equation,  Equation  (6)  is  given 
as 


a  =aT  +Kcos  (cdt+8) 


(16) 
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where,  for  quadratic  variations  in  Equation  (15). 


k=2ko{cm^cm^/1(cm^cm^  mz 

+4CCm^CM4)o]}2 


Equation  (16)  may  be  fit  to  the  angular  data  and  Equations  (17)  to  (18) 
then  employed  to  obtain  the  non-linear  variation  of  the  stability  coefficients 
with  angle  of  attack.  Figure  8  shows  C*w,'«*)  to  be  non-linear  and  of  a’feoft- 
spring"  nature.  Figure  9  shows  that  +  Cm*  also  demonstrates  non¬ 

linearity,  but  its  variation  is  of  a  "hard  spring’nature . 

It  has  been  shown3,4,8  that  for  the  case  of  a  variable  frequency 
(Crt«/  0  in  the  wind  tunnel)  the  amplitude,  K(t),  must  be  corrected  as 
follows 

WKCO^F  <20> 

In  order  to  examine  the  degree  of  accuracy  of  the  approximate  non-linear 
solution,  numerical  integrations  of  the  equation  of  motion,  Equation  (6), 
were  performed  using  a  Fourth  Order  Runge  Kutta  scheme.  Nominal  values 
of  the  aerodynamic  coefficients  were  input  and  the  subsequent  time  history 
of  «*(t)  was  then  fit9  and  the  non-linear  theory  applied  to  the  results.  Table 
II  shows  the  percent  errors  obtained  from  this  analysis. 
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TABLE  II.  PERCENT  ERROR  FROM  ANALYSIS 


Coefficient 

Input 

Output 

Percent  1 

Lm*„ 

-30.0 

-29.94 

0.18 

+  15.0 

+15.07 

0.47 

-75.0 

-74.83 

0.21 

Cm.j'*'  C r.i,  )* 

-50.0 

-55.74 

11.49 

Cm  % 

1.0 

1.0 

0. 

Vk 


A  model  of  the  basic  finner  was  mounted  on  a  jewel  support  to 
minimize  bearing  friction  in  the  angular  oscillations.  ?  Again,  optical  techni¬ 
ques  were  employed  to  determine  the  angular  data.  Several  tests  were  con¬ 
ducted  at  one  c.g.  position  to  establish  repeatability  of  the  stability  coeffi¬ 
cients. 

This  dramatic  difference  in  the  nature  of  the  non-linearities  of  these 
two  important  stability  coefficients  requires  further  study. 

Of  particular  interest  is  the  high  degree  of  repeatability  of  the  test 
runs.  This  fact  provided  added  confidence  in  the  application  of  the  non-linear 
theory. 


Fig.  8.  Cm*  versus  <x 
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Fig.  9.  CMq  +  Cj^  versus  a 


Mean  values  of  the  coefficients  of  the  polynomial  expression  for  the 
aerodynamic  stability  coefficients  of  Equation  (15)  are  shown  below. 


^  M 


=  -  22. 65 


=  68.31 


CM,+CHfc).  =  -  193.58 


(C,,  +  Cm,A=  -  19334. 


These  results  suggest  that  the  non-linear  method  of  analysis  should  be 
used  in  order  to  obtain  accurate  values  for  the  stability  coefficients.  It  has 
been  shown4  that  the  "log  -  decrement"  technique  yields  highly  inaccurate 
results  for  the  damping  moment  when  the  restoring  moment  is  non-linear. 
Although  use  of  Equation  (2u)  tends  to  provide  better  results  in  this  regard, 
it  is  strongly  recommended  that  the  complete  non-linear  method  be 
employed. 


btSUlIO'N  IV 


ROLL  LOCK-INv 


It  is  well  known  that  missiles  and  sounding  rockets  experiencing 
resonance  instability  have  suffered  abnormally  large  angles  of  attack  un¬ 
explainable  by  resonance  alone.  This  phenomena,  known  as  Catastrophic 
Yaw  ,  results  from  a  lunar -type  motion  of  the  vehicle  in  the  resonance 
region.  This  motion  is  characterized  by  the  angle  between  a  given  fin  and 
the  (rotating)  plane  of  the  angle  of  attack  being  essentially  constant,  i.e. 
roU  lock-in.  3  in  order  to  more  fully  understand  the  effects  and  hopefully 
control  roll  lock-in,  the  simplified  case  of  pure  rolling  motion  at  a  constant 
angle  of  attack  was  first  examined . 

The  governing  equation  for  a  rolling  missile  at  constant  angle  of 
attack  i « 

IxP  =  L(Sa)  +  LCp^a)*  LC^a) 


where 


L(  S  a.)  =  CA(  8  a  )  Q  S  d  *  ro11  moment  due  to  fin  cant 
L(p,y,a)  =  C/  P, y,a)  Q  Sd  =  ro11  dampin«  mornent 

L(y,a)  =  Cj^cOQSd  *  induced  roll  moment 

p  =  y  for  a  =  constant. 

Theory  and  experiment  have  shown  that  the  induced  roll  moment  for  a 
cruciform  finned  missile  may  be  approximated  by  a  sine  function  as 

L( y}a)  =  L(a)  sin  4y 

Therefore  the  differential  equation  of  motion  may  be  written 

LCSa)  *- LCp,7,a)  +  L(a)  si  n  4y  =  a 

Linear  Analysis 

For  the  case  where  the  roll  damping  moment  is  assumed  to  vary 
linearly  with  the  rolling  velocity  and  the  inherent  non-linear  sine  function 
approximated  by  a  straight  line  tangent  to  the  curve  at  the  trim  angle, 
Equation  (21)  may  be  written 


{ 

5 


(22) 


LjaSa  +  LpC a)y  +  Lyja  =  Ixy 


Since  the  coefficients  of  the  differential  equation  are  assumed  constant,  the 
solution  for  Equation  (22)  is  the  same  as  the  solution  for  pitching  motion  dis 
cussed  previously,  and  therefore  we  may  write 


2IxX 

*P  ^QSd 


(23) 

(24) 

(25) 


Two  horizontal  wind  tunnel  tests  were  conducted  at  supersonic 
speeds  10  to  obtain  a  time  history  of  'Y  v  for  a  Basic  Finner  exhibiting  roll 
lock-in  at  a  constant  angle  of  attack.  One  test  demonstrated  damped  roll 
oscillations  while  the  other  undamped.  Figures  10  and  11  show  the  damped 
and  the  undamped  oscillation  respectively. 


TlMGCscc-l 

Figure  10.  Roll  Orientation  versus  Time 
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Figure  11.  Roll  Orientation  versus  Time 


By  fitting  Equation  (9)  to  the  above  angular  data  by  the  Method  of  Differential 
Corrections  to  obtain  \  and  GJ  ,  Equations  (24)  and  (25)  were  used  to  com¬ 
pute  Cjp  and  Cgr  .  These  results  are  presented  in  Table  III. 


TABLE  III.  STABLE  AND  UNSTABLE  ROLL  COEFFICIENTS 

/N/ 

Dynamic  Stability  Mach  No.  Co  Ct 

Stable  3.50  -6.66  -.027 

Unstable  3.50  1.11  -.019 


Non-Linear  Analysis** 

The  first  approximation  of  Kryloff  and  Bogoliuboff 12  was  applied  to 
the  non-linear  differential  equation  of  motion,  Equation  (21), for  the  case  of 
zero  fin  cant.  While  the  induced  roll  moment  term  inherently  non-linear, 
the  functional  dependence  of  the  roll  damping  moment,  Up,  on 
was  also  included. 

By  considering  qualitatively  the  damping  force  distribution  across  the 
fins  combined  with  the  fin-body  wake  effect,  the  dependence  of  L  ( p.  *,'*) 
on  may  be  determined.  As  in  the  case  of  the  linear  analysis  L(p,s  ,cx) 
was  assumed  to  vary  linearly  with  the  rolling  velocity. 


I 

! 
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L(p,y,a)  -Lp(-y,a)p  -  (Lp^a)+Lpfe)(y-  ^)2}p<26) 


7 


Therefore  the  equation  of  motion  may  be  written 

-j-^-SIN  4y  =  0.(27) 


I 


X  *x  J>x 

An  approximate  solution  to  the  above  non-linear  equation  is  given  as 


7  =  K  C0S(wt  +  8) 


(28) 


where 


K  =  2K(0)[C  /(4C  ♦ 

r  O  “o 

*  2  r  CoQSd*^  n  j/2  Co  QSd1 

QpK(0)  {l-e  21. V  tJ)J  g  *41 \y  t 
w  =  ?  +  c  (oIT  -^.4K)- 

I*  ‘  v  I*  *  K 


(29) 


(30) 
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In  addition,  the  non-linear  damping  moment  enabled  the  accurate  prediction 
of  a  limit  cycle  in  roll  oscillation.  The  non-linear  solution  yields  an  equa¬ 
tion  for  the  time  rate  of  change  of  the  amplitude  as  follows 


dK  LP.  ,,  .  LPl  3 

dt  2 Ix  *  8IX  * 


Hence ,  for  a  limit  cycle 

dK 

dt 


and  therefore 


o 


2V~Q,,./c 


=  constant 


(31) 


LP„<  0  lp>  0 

To  investigate  the  accuracy  of  the  limit  cycle  prediction,  Equation  (27)  was 
numerically  integrated  with  the  linear  and  non-linear  damping  moments  re¬ 
lated  to  the  initial  amplitude  as  given  by  Equation  (31).  Results  yielded  roll 
oscillations  which  damped  0.25°  after  10  cycles  of  oscillation. 


In  addition,  numerical  studies  were  performed  to  substantiate  the 
approximate  non-linear  solutions,  Equations  (28)  to  (30).  The  governing 
differential  equation  was  numerically  integrated  to  yield  a  time  history  of 
*Y  It)  .  Equation  (28)  was  then  fit  to  these  results  by  the  Method  of  Differen¬ 
tial  Corrections  and  Equations  (29)  and  (30)  applied  to  compute  C  *  Po  , 

C  ,  Cjj  (<*)  .  The  percent  errors  were  computed  with  respect  to  the 

input  values.  A  correction  to  the  amplitude,  similar  to  that  of  the  1-D 
pitching  motion,  must  be  applied  due  to  the  frequency  variation.  The  results 
shown  are  for  a  damped  roll  lock-in  oscillation  with  an  initial  roll  angle  of 
70°.  See  Table  IV. 


TABLE  IV.  PERCENT  ERROR  OF  COEFFICIENT  COMPUTATION 


Coefficient 

Input 

Output 

%  Error 

CXp, 

-2.0 

-2.02 

1.00 

C*p 

15.0 

14.42 

3.87 

C,  M 

0.5 

.4998 

.04 

Subsonic  wind  tunnel  tests  are  currently  being  undertaken  to  obtain 
aerodynamic  data  on  the  Basic  Finner  configuration  from  angular  measure¬ 
ment  of  roll  oscillations.  A  jewel  support  system  mounted  in  a  vertical 
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wind  tunnel  at  the  Aero-Space  Engineering  Department  Univer- 
sity  of  i'iotre  Dame,  has  drastically  reduced  the  effect  of  bearing  friction. 

At  the  present  time  calculations  have  been  made  to  determine  the  induced 

coefficient  at  several  angles  of  attack.  The  results  shown  below 
are  for  a  basic  finner  model  with  flat  plate  fins. 


0-6 


>5 
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Basic  Fiw.mer 
(PLAT  PLATE  Flin) 


o 

Figure  13. 
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The  variation  of  the  induced  roll  moment  coefficient  with  a  com¬ 
pares  feivorably  with  that  of  Rhodes  and  Shannon  13  on  cruciform -finned  con- 

T*?  ’  In  a^dltion’  the  data  in  FiSure  demonstrates  the  possibility 
that  CjW  may  be  negative  for  a<  ,  .  This  is  also  in  agreement  with  die 

results  obtained  on  other  cruciform  finned  configurations.  1'* 


■» 
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SECTION  V 


THREE  DEGREE  OF  FREEDOM  MOTION 


A  three  degree  of  freedom  pitching,  yawing,  and  rolling  model  was 
designed  for  low  subsonic  dynamic  wind  tunnel  testing.  These  tests  included 
both  A  ball  bearing  side  support  system  for  use  in  a  horizontal  tunnel  14  and 
a  jewel  cup  rear  support  system  for  use  in  a  vertical  "blow  down"  tunnel. 15 
In  the  case  of  the  side  support,  the  entire  model  was  free  to  roll  except  for 
a  small  cylindrical  section  at  the  center  of  gravity,  while  for  the  jewel 
support  system  no  rolling  constraints  were  applied. 


The  differential  equation  for  3-D  angular  motion  may  be  written  as 


Maa  +M  t  +  M^a  +  i  MpaPa  =  15  -  i  PI  a  (32) 


The  general  solution  when  the  coefficients  are  constant  may  be  written  as 


a  =  K  ec +  k  eiPt 


where 


\,-(s 


Co  ~  pt 
1,2  Hi 
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^CMd}  r  |  +T'i  +  ^MPar  I  QSd 
I  U~  J  ‘  Ix  ;  2V 


t  =7TTT^ 

0  +  ^  4ICMaQSd 

The  aerodynamic  coefficients  may  then  be  solved  as 


0  + 


c  = 

**  Q  Sd 


(36) 


? 


S^+^Md 


'Mpa 


(  \  07  v  / 

v.'V*2-'  CIM 

(37) 

QSd2 

• 

^  ^  +  ^Md 

(38) 

Q  Scpr  21 

These  tests  were  directed  towards  an  evaluation  of  the  stability 
oefficients  from  3-D  angular  motions  with  primary  emphasis  on  obtaining 


r 


In  a  similar  manner  as  for  the  1-D  pitching  motion,  optical  techniques 
were  employed  to  obtain  the  experimental  3-D  angular  data.  An  example  of 
the  complex  angular  motion  is  presented  in  Figure  14. 


By  overlapping  fits  of  small  sections  of  data  the  motion  parameters  can  be 
determined  as  continuous  functions  of  time .  The  aerodynamic  stability 
coefficients  may  then  be  determined  from  these  motion  parameters. 
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Sectional  fits  of  Equation  (33)  to  the  above  data  resulted  in  a  "probable 
i  error  of  fit"  ranging  from  0.5°  to  0.2°  as  shown  in  Figure  15.  The  resulting 
mean  percent  error  is  approximately  2%  to  3%. 
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Figure  15.  Probable  Errors  versus  Time 

The  motion  parameters  obtained  from  sectional  fits  of  the  angular 
motion  of  Figure  14  are  presented  as  functions  of  time  in  Figures  16  to  18. 

It  is  seen  from  these  results  that  small  non-linear  variations  are 
present. 

Using  these  motion  parameters  in  addition  to  wind  tunnel  conditions 
and  model  mass  parameters  the  stability  coefficients  Cm*  ,  CMs  +  Cm* 
Cmp*  were  computed  from  Equations  (36),  (37),  and  (38)  respectively. 
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Figure  16.  Kj  and  K2  versus  Time 
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Figure  17.  Frequencies  versus  Time 
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The  test  runs  were  conducted  at  low  subsonic  speeds  at  one  c.g.  position  . 


Figures  19,20,21  present  (V  .  C 
mean  angle  of  attack.  The  sectional 


m*  'r  >  ''tip*  as  a  function  of 

fits  yielded  variations  of  the  stability 


coefficients  with  angle  of  attack  which  compare  favorably  with  those  obtained 
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Figure  19.  Cm  versus  a 
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Figure  20.  C^q +  PMq  versus  a 
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Fig.  21 
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Figure  21.  Cmd  .  versus  a 
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In  addition,  o-u  tests  were  umuuticu  w..  -  *  ***“ — -  - 

on  a  jewel  sting  support  in  i,he  Notre  Dame  vertical  "blow  down  subsonic 

wind  tunnel.  The  3-D  oscillations  were  obtained  and  reduced  in  a  manner 

similar  to  that  of  the  above  analysis.  It  should  be  noted  that  this  is  a  rear 

support,  in  contrast  to  the  side  support  of  the  ball  bearing  system.  This 

jewel  system  minimizes  the  bearing  friction  as  welj  as  reducing  the  support 

interference.  The  stability  coefficients  CM,  and  CMp(?  obtained  from 

several  tests  are  presented  in  Figures  22  and  23.  For  the  tests  conducted, 

only  the  precession  mode  of  oscillation  was  obtained.  Hence,  a  value  of 

r  •+  r  =  -  350  was  used  in  determining  the  values  of  CM 

v*  m  -  Vt  i  r 
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Figure  22.  CM  versus  a 
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Figure  23.  C^p  versus  a2 


In  the  comparison  of  C  M  A  it  should  be  noted  that  the  values  in  Figure 
22  are  larger  than  those  of  the  previous  1-D  and  3-D  analyses.  This  results 
from  the  c.g.  location  in  addition  to  the  fact  that  the  jewel  supported  model 
was  equipped  with  roll  tabs  on  the  fins ,  thereby  increasing  the  fin  area  by 
12.5%. 


The  various  stability  coefficients  of  this  linear  analysis  appear  to  vary 
with  angle  of  attack.  Thus  it  is  essential  to  reconsider  the  solution  of  the 
equation  of  motion  in  which  the  stability  coefficients  are  non-linear  functions 
of  the  angle  of  attack. 


Non-Linear  Analysis 

For  the  case  where  the  stability  coefficients  are  non-linear  functions 


of  angle  of  attack 
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(39) 


„,2 


C„  Cl*l)  =  C„  +  Cm  I  cT  I  + 

^pa  ^pa,  Mpc^1  1 

an  approximate  solution  for  the  complex  angle  of  attack  is  given  as 


a  =/?  +ja  =  K,e 


CX,+  iw.)t 


+  K2e' 


CX5+ iaOT 


where  for  quadratic  variations  in  Equation  (39), 

>  .  QSd2/2IVC«a2  ,  '  2  , 

72~  2V  l^cpg^ 
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was  conducted.  Sample  aerodynamic  stability  coefficients  in  the  form  of  the 
quadratic  variations  in  Equation  (39)  were  input  to  a  computer  program  which 
numerically  integrates  the  6-D  equations  of  motion.  A  wind  tunnel  angular 
motion  was  simulated.  Using  the  a  and  output  by  the  6-D  computer 
program,  fits  of  the  theory,  Equation  (40),  were  made  to  the  angular  data  using 
the  WOBBLE  computer  program  to  determine  the  stability  coefficients.  A 
comparison  of  the  results  obtained  from  the  WOBBLE  fits  with  those  input  to 
the  6-D  program  are  given  in  Table  V.  As  shown  in  these  results,  the  approx¬ 
imate  non-linear  solution  accurately  represents  the  numerically  integrated 
data. 

This  technique  was  applied  to  several  tests  of  a  particular  3-D  of  free¬ 
dom  cruciform-finned  model  ai.  ’ow  subsonic  velocities.  "Fits"  of  Equation 
(40)  through  (42)  yields  the  non-linear  aerodynamic  stability  coefficients. 
Figures  24  through  26  present  the  non-linear  variation  of  CM(1  ,  Cn  +-  CMa 
and  CM  with  angle  of  attack,  01  * 


TABLE  V. 

EVALUATION  OF  NON-LINEAR  SOLUTION 

6-D 

Wobble 

Coef. 

Input 

Output 

%  Error 

CmC1o 

-113. 

-113. 

0 

r 

1067. 

1060. 

0.7 

VMt  LMc 

j  "300°* 

-3006. 

0.2 

)  49249. 

50570. 

2.7 

/  . 

^  M 

1000. 

1003. 

0.3 

16416. 

15780. 

3.9 

i 
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Figure  24.  versus  a 


Mean  values  of  the  coefficients  for  the  polynomial  representation  of 
the  angle  of  attack  variation  of  the  stability  coefficients  were  found  to  be 


1122. 

212,633. 

6469. 
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In  analyzing  the  above  results  it  is  seen  that  the  restoring  moment 
exhibits  "soft -spring"  characteristics  while,  the  damping  moment  shows  a 
"hard  spring"  variation.  This  is  the  same  trend  as  observed  in  the  previous 
1-D  and  3-D  results.  Further  study  is  required  in  this  regard. 


SECTION  VI 


CONCLUSIONS 


Various  non-linear  methods  of  data  analysis  have  been  developed  and 
evaluated  by  utilizing  computer  generated  data  and  by  comparing  the  results 
with  the  original  inputs.  In  the  case  of  the  non-linear  pure  pitching  motion, 
C„„  and  (  C*H  +  C*A  )  agree  to  .  18%  and  .21%  respectfully,  and  the  non¬ 
linear  quadratic  coefficients  and  (C^q  +  CM<i  )2  agree  to  .47%  and 
11.5%.  In  the  case  of  the  nonlinear  pure  rolling  motion  Cx?o  agrees  to  an 
accuracy  of  1.0%  and  the  non-linear  quadratic  coefficient  C/p?  (u*)  and 
C  n  M  agree  to  3. 87%  and  .  04%  respectfully.  In  the  case  o?f  combined 
pitching,  yawing  and  rolling  motion  and  (  +  C^)  agree  to 

.01%  ana  .2%.  The  hbn-linear  coefficients  C*„t  and , c\.  +  C«J,a8ree  t0  •?% 

and  2%.  Of  particular  importance  is  the  Magnus  coefficient  C  which 
agrees  to  .3%  and  the  non-linear  quadratic  Magnus  coefficient,  *peUCft. 
which  agrees  to  3.9%. 

Various  new  dynamic  wind  tunnel  testing  techniques  were  developed  and 
the  data  was  reduced  by  using  the  new  non-linear  methods  of  analysis.  The 
experimental  data  was  found  to  be  well  represented  by  the  theories  for  the 
three  non-linear  motions  considered:  pure  pitching,  pure  rolling  and  com¬ 
bined  pitching,  yawing,  and  rolling.  The  values  for  the  associated  stability 
coefficients  *  <  C*,  +  C**  )»  and  C*  («)  compare  well 

in  repeated  tests  and  the  detailed  results  are  all  represented.  It  is  concluded 
that  the  various  motions  are  definitely  non-linear  and  that  the  non-linear 
methods  must  be  used  for  their  reduction  and  analysis.  It  is  suggested  that 
previous  disagreement  between  data  from  various  facilities  may  be  due  to 
these  serious  non-linear  effects. 

For  the  Basic  Finner  configuration  it  i?  found  that  the  non-linear  re¬ 
storing  moment  is  of  the  "soft -spring"  type  while  the  non-linear  damping 
and  lag  moments  are  of  the  "hard  spring"  type.  Future  study  is  suggested  in 
this  area. 

Of  particular  importance  are  the  new  dynamic  wind  tunnel  testing 
techniques  utilizing  special  jewel  supports  which  permit  nearly  frictionless 
free  motion.  The  unique  vertical  "downdraft"  wind  tunnel  has  yielded  excel¬ 
lent  data  on  pure  rolling  motion  and  on  combined  pitching,  yawing,  and 
rolling  motions . 


33 

(The  reverse  of  this  page  is  blank) 


REFERENCES 


1.  Bolz,  R.E.,  and  Nicolaides,  J.D. ,  A  Method  of  Determining  Some  Aero¬ 
dynamic  Coefficients  from  Supersonic  Free -Flight  'tests  of  a  Rolling 
Missile,  Journal  of  the  Aeronautical  Sciences,  Vol.  17,  No.  10,  6ct. , 
1950,  BRL  Rept.  No.  711. 

2.  Nicolaides,  J.D. ,  On  the  Free  Flight  Motion  of  Missiles  Having  Slight 
Configurational  Asymmetries,  BRL  ftept.  No.  856,  1953,  and  IAS 
Preprint  Isto.  395,  1952’. 

3.  Nicolaides,  John  D. ,  Free  Flight  Dynamics ,  Text,  Aero-Space 
Engineering  Department,  University  of  Notre  Dame,  1961. 

4.  Eikenberry,  R. ,  Wobble,  Analysis  of  Missile  Dynamic  Data,  Part  I: 
Angular  Motion,  Aero-Space  Engineering  Department,  University  of 
Notre  Dame,  1967. 

5.  Nelson,  R.C. ,  On  the  Pure  Pitching  Motion  of  the  Basic  Finner  Missile, 
Master's  Thesis,  tJepartment  of  Aero-Space  Engineering,  University  of 
Notre  Dame,  1966. 

6.  Nicolaides,  John  D. ,  Eikenberry,  R.S.,  Dynamic  Wind  Tunnel  Testing 
Techniques ,  AIAA  Aerodynamic  Testing  Conference,  Los  Angeles, 
California ,  1966. 

7.  Clare,  T.A.,  Chin,  A.B. ,  Lando,  D.W.,  A  Technique  for  Obtaining 
Non-Linear  Aerodynamic  Stability  Coefficients  for  a  Missile  Pure  Yaw, 
Aero-Space  Engineering  Apartment  Report,  University  of  Notre  Dame, 
1968. 

8.  Murphy,  C. ,  Free  Flight  Motion  of  Symmetric  Missiles ,  BRL  Report 
No.  1216,  Aberdeen  Proving  Ground,  Maryland,  1963. 

9.  Clare,  T.A. ,  FIT,  A  Program  to  Fit  Analytic  Functions  to  Data  by 
Differential  Corrections,  Aero-Space  Engineering  Department  Report, 
University  of  Notre  Dame ,  1968 . 

10.  Bozzonetti,  E.,  Analysis  of  the  Oscillating  Motions  Associated  with  Roll 
Lock-In,  Master's  Thesis,  Department  of  Aero-Space  Engineering, 
University  of  Notre  Dame  ,  1966. 

11.  Clare,  T.A. ,  A  Non-Linear  Analysis  of  Dynamic  Roll  Lock-In  in  One 
Degree  of  Freedom,  Aero-Space  Engineering  Department  Report, 
University  of  Notre  Dame,  1968. 


35 


12.  Klyloff.  N. .  Bogolinhoff,  N. ,  Introduction  to  Non-Linear  Mechanics, 
Translated  by  S.  Lefschetz,  Princeton  Univ.  Press,  Princeton,  New 
Jersey,  1947. 

13.  Rhodes,  C.W. ,  Shannon,  J.H.  W. ,  Results  and  Conclusions  of  the 
Joint  R.  A .  E .  /W .  R .  E .  Research  Programme  on  the  Flight  bynamics 
and  Ballistic  Consistency  of  Freely  Falling  Missiles  ,  Part  i:  kombs 

,  Stabilized  by  Fixed  Cruciform  Fins,  Department  of  S upply , A ustralian 
Defense  Scientific  Service,  Weapons  Research  Establishment,  HSA 
20,  1965. 

14.  Lorenzen,  C. ,  On  the  Determination  of  the  Magnus  Moment  Coefficient 
for  the  Basic  Finner,  Master’s  Thesis,  Department  of  Aero-Space 
engineering',  University  of  Notre  Darne,  1967. 

15.  Foust,  J. ,  On  the  Investigation  of  a  Vertical  Wind  Tunnel  System  for 
Determining  Aerodynamic  Moment  Coefficients,  Masters  Thesis, 
Department  of  Aero-Space  Engineering,  Univer s ity  of  Notre  Dame , 
1968. 

16.  Ingram,  C. ,  An  Approximate  Solution  of  the  Non-Linear  Differential 
Equation  for  the  Complex  Angle  of  Attack  of  a  Symmetric  Missile  , 

Air  Proving  Ground  Center,  Eglin  Air  Force  Base,  Florida.  (Pending 
Publication). 

17.  Ingram,  C. ,  Martin,  J. ,  A  Three  Degree  of  Freedom  Subsonic  Wind 
Tunnel  Test  of  the  Apache  Sounding  Rocket,  White  Sands  Missile  Range, 
New  Mexico  (Pending  Publication). 

18.  Nicolaides,  Eikenberry,  Ingram  and  Clare,  On  the  Flight  Dynamics 
of  the  Basic  Finner  in  Various  Degrees  of  Freedom.  AIAA  Paper 
No.  68-890. 


t 


36 


UNCLASSIFIED 

DOCUMENT  CONTROL  DATA  *  R  &  P 

_ _ (Security  cImbm  Iflcmtion  ot  tltie,  body  of  abstract  end  Indexing  annotation  mutt  fee  entered  whe  n  fb#  o  vtnll  repott  fo  c  tereltted^ 

1  ORIGINATING  ACTIVITY  ('Co/poraf#  •ufhor;  2a*  REPORT  SECURITY  CLASSIFICATION 

Department  of  Aero-Space  Engineering  UNCLASSIFIED 

University  of  Notre  Dame  **>•  sroup 

Notre  Dame ,  Indiana _ _ 

FLIGHT  DYNAMICS  OF  THE  BASIC  FINNER  IN  VARIOUS  DEGREES  OF 
FREEDOM 

4,  DESCRIPTIVE  NOTES  (Typ»  ot  report  end  inclueive  deter) 

Final  Report  -  February  1966  to  May  1968 _ 

5.  AUTHORISE  (Fitef  nemet  middl e  tnttiet,  teet  nmme)  ~ 

John  D.  Nicolaides  Charles  W.  Ingram 

Roberts.  Eikenberry  Thomas  A.  Clare 


|0.  REPORT  DATE 


1 7a.  TOTAL  NO.  OF  PAGES 


76.  NO.  OF  REFS 


July  1968 _ 

•a.  CONTRACT  OR  GRANT  NO. 

AF  08  (635)- 5275 

b.  PROJECT  NO. 


Ida.  ORIGINATOR'S  REPORT  NUMSER(S) 


9b.  OTHER  REPORT  KO(S)  (Any  other  number*  thmt  tnmy  be  meetjnod 
thle  report) 


AFATL-TR-68-82 


10.  DISTRIBUTION  STATEMENT 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of 
the  Air  Force  Armament  Laboratory  (ATBR),Eglin  AFB,  Florida  32542. _ 


|  1 1  •  SUPPLEMENTARY  NOTES 


j  i  2  ■  SPONSORING  MILITARY  ACTIVITY 


Available  in  DDC 


Air  Force  Armament  Laboratory 
Air  Force  Systems  Command 
Eglin  Air  Force  Base,  Florida 


113.  ABSTRACT 


Basic  Finner  models  are  tested  in  both  horizontal  and  vertical  wind  tunnels  in  order 
to  develop  dynamic  testing  techniques  and  to  perfect  new  metnodb  of  non-linear 
analysis.  The  simple  free  motions  of  pure  pitching  and  pure  rolling  are  first  studied. 
Then,  free  angular  motions  of  combined  pitching,  yawing,  and  rolling  are  investi¬ 
gated.  While  the  linear  analysis  methods  appear  to  yield  good  results,  the  new  non¬ 
linear  analysis  methods  are  found  to  be  essential  in  providing  an  accurate  repre- 
:  sentation  of  the  various  motions  and  in  providing  correct  values  for  the  various 
j  static  and  dynamic  stability  coefficients ,  ^  > 


FORM 

\  NOV  68 


1473 


UNCLASSIFIED 


Security  Classification 


UNCLASSIFIED 


Security  Classification 


rr  ■  - - - - 

I  KEY  WORD* 

LiNn  A 

LINK  a 

LINK  C 

ROLE 

WT 

(tout  I  WT 

1 

Flight  Dynamics 

Basic  Finner 

Degrees  of  Freedom 

Non*  Linear  Analysis 

Linear  Analysis 

UNCLASSIFIED 


Security  Classification 


